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The question how to introduce thermal fluctuations in the equation of motion of a magnetic 
system is addressed. Using the approach of the fluctuation-dissipation theorem we calculate the 
properties of the noise for both, the fluctuating field and fluctuating torque (force) representation. 
In contrast to earlier calculations we consider the general case of a system of interacting magnetic 
moments without the assumption of axial symmetry. We show that the interactions do not result 
in any correlations of thermal fluctuations in the field representation and that the same widely used 
formula can be used in the most general case. We further prove that close to the equilibrium where 
the fluctuation-dissipation theorem is valid, both, field and torque (force) representations coincide, 
being different far away from it. 



The problem of a correct introduction of temperature 
in the equation of motion of a magnetic system has gained 
much importance as a result of technological require- 
ments of magnetic recording industry ||, ||, [|. This is 
associated with the need to perform calculations of mag- 
netization dynamics at finite temperatures. Open prob- 
lems include fast magnetization switching, thermal sta- 
bility and magnetic viscosity, among others. The correct 
solution of the problem is still far from being understood. 
The main difference between the magnetic problem and 
the standard molecular dynamics approach is that the 
magnetic moment dynamics is governed by the Landau- 
Lifshitz equation which includes the precession of a mag- 
netic moment around its internal field direction. It com- 
prises coupled first-order equations for the magnetization 
components, and the requirement of conservation of the 
magnetization amplitude. As a consequence, no analogue 
of mass and kinetic energy exist in the system, thus mak- 
ing it impossible to introduce the temperature through 
this mechanism. 

Consequently, the temperature is introduced through 
small deviations from the equilibrium configuration. 
Therefore, strictly speaking, this approach is only valid 
when these deviations are small and it cannot be used 
for fast magnetization switching. 

Let us briefly summarize the original approach from W. 
Brown [Q, The underlying equation of motion is the 
Landau-Lifshitz-Gilbert equation which can be written 
in the form 
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70 is the gyromagnetic ratio and a is the damping con- 
stant. The magnetic moment rvJ is normalized to the 
saturation value M s , and the internal field S is normal- 
ized to the anisotropy field Hk = 2K/M S . The energy 
E* = E/2KV, where K is the anisotropy value and V 
is the particle volume, contains all the necessary energy 
contributions: anisotropy, exchange, magnetostatic and 
Zeeman. 

W. Brown proposed the inclusion of thermal fluctua- 
tions via a random field, added to the internal field, Eq. g. 
For the calculation of the properties of the random field 
he outlined two methods: (i) based on the fluctuation- 
dissipation theorem (see also || ) and (ii) by imposing the 
condition that the equilibrium solution of the correspon- 
dent Fokker-Planck equation is the Boltzmann distribu- 
tion (see also ) . As a result of both the thermal field 
statistical properties are given by 
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where i,j denote Cartesian components x, y, z. Different 
approaches based, for example, on the Landau-Lifshitz 
rather than on the Landau-Lifshitz-Gilbert equation were 
also introduced M |). 

However, the properties of the thermal noise, Eqs. |^, 
were derived only for one isolated particle. Moreover, 
Brown considered in his paper only the simplest axi- 
ally symmetric case. Nevertheless, in the past the formu- 
las above provided the basis for practically every numer- 
ical method ||, 0, [l0| |llj for the computation of mag- 
netization dynamics taking into account thermal fluctu- 
ations. But the investigated magnetic systems usually 
comprise interacting particles |l2|, |l3|, |l4|, Kjk |ll| [l7| due 
to magnetostatic and/or exchange couplings. For that 
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case the thermal field may be expected to be influenced 
by correlations between different particles ]l0| . Hence, 
it is necessary to generalize Brown's result to the case of 
interacting magnetic moments. To the best of our knowl- 
edge, this has never been done before. 

In what follows we start with the Brownian dynamics 
approach (see JTa]) which was originally applied to mag- 
netic systems by A. Lyberatos et al. || [H^. However, we 
consider the general case of an interacting system with a 
non-axially-symmetric potential. Following the standard 
approach, we introduce the temperature into the motion 
of the Brownian particles (i. e. the magnetic moments) 
as a result of the fluctuation-dissipation theorem. Conse- 
quently, this approach is only valid when small deviations 
from equilibrium are considered. 

The general Langevin equation of motion is written in 
the form 



dS 



(4) 



where the jij are the so-called kinetic coefficients, Xj are 
variables which are thermodynamically conjugate to Xj, 
and S is the entropy of the magnetic system. For a closed 
system in an external medium, 



x _ 1 d_E 

3 k B T dx. 
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In Eq. |^, fi is a random force representing thermal 
fluctuations in the system having the properties 



(/*(*)) = and (fi(0)f j (t))=t Hj 6(t) 
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where 



A linear equation of motion of the form 
dx 
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with the associated energy 
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can be rewritten as 

dxj 
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so that the matrix is related to the kinetic coefficients 
7ife in the following way H : 



L ik - "j-yE' 1 '^ 



(11) 



In micromagnetics the motion of a magnetic moment 
M is governed by the deterministic LLG equation (Eq. [I]) . 
For the equilibrium state of the system Brown's condition 



M? x Hf = 



(12) 



must be satisfied, implying that here and M. i are 
parallel. Close to equilibrium, the LLG equation can be 
linearized using small deviations 



= Mi - m' 



H, 



from their equilibrium values, yielding 

, 37V 

drrii x - 



(13) 
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Here, the indices i, j count the particles sites 1,...,N as 
well as their x,y,z coordinates. The internal fields hj 
play the role of the variables which are thermodynami- 
cally conjugate to mj, 



1 dE 2KV h 
3 ~ t B 1 Td^n~ ~ ~ k B T j ' 



(15) 



Thus, the LLG equation should be rewritten in the 
form 



drrii 
~df 
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which is an easier way to calculate the kinetic coefficients 
than the use of Eq. |ll|. The representation of the LLG 
equation in the form of Eq. ^] means that in what follows 
the thermal fluctuations are introduced as a fluctuating 
torque (a generalized force rather than a field). Later we 
will show that in the linear approximation this is equiva- 
lent to the standard fluctuating field representation. Al- 
ternatively, Eq. |l6| could be viewed as a polar represen- 
tation of the magnetization vector m\ = 0i,rnf = ipi, in 
this case the conjugate variables are the polar projections 
of the internal fields (he, h v ) and the fluctuations fi will 
stand for the random field polar components. This latter 
approach was used originally by W. Brown [4j. 

We continue by writing the energy of the system in the 
form 



N 
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where A is the Lagrange multiplier. In the zero order 
approximation one obtains 



1 X 1 



(18) 



which corresponds to Brown's condition, Eq. The 
linear approximation leads to the equilibrium condition 



M i'- ii'; 
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which leaves us only the quadratic form for the energy 
expression near the equilibrium, 



E* = E 



N 



2 l} 



(20) 



The general expression for magnetic energies is a 
quadratic form in terms of the magnetization (apart from 
the Zeeman term which is included in the equilibrium 
field Ilf and condition |l^) . Therefore, it is reasonable to 
suppose that the total field can be expressed as 



= h 



eff,i 



A? 



(21) 



where ft."^ i are the components of the effective field due 
to the different energy contributions and — Xmf is the 
field due to the kinematic interaction expressing the con- 
straints (Lagrange multiplier) . The value of the Lagrange 
multiplier is normally found from the equilibrium condi- 
However, its actual value is not necessary for 
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tion 

calculations due to the fact that the LLG equation con- 
serves the magnetization length. Latin indices represent 
the sites of the moments and the Greek ones the magne- 
tization components x, y, z. In this case the final expres- 
sion for the energy (Eq. EG) takes the form 



E* = E - ]T (B 



a/3 
ij 



>\j>v, ,;•'»'/'»;. (22) 



The expressions for the kinetic coefficients could be 
obtained by using directly the expression [ll]. In this ap- 
proach it seems that the final result could also include 
correlations between different particles |fl0|| . But this 
is not the case: the kinetic coefficients can be obtained 
much easier representing the linearized LLG equation in 
the form of Eq. [l^ yielding 
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(24) 
(25) 
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The other coefficients can be obtained by symmetry. 
Note, that there are no correlations between different 
particles. Also, the kinetic coefficients have obviously 
reversible parts (coming from rotation) and irreversible 
parts (from damping). The reversible antisymmetric 
parts do not contribute to the thermal fluctuations af- 
ter adding the kinetic coefficients to calculate the matrix 
\i from Eq. 0, yielding 

ak B T [, 



t'i 
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(27) 
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Once again, the others can be obtained by symmetry. 
Note that in a general system of coordinates there are 
correlations between different magnetization components 
but no correlations between different particles. However, 
if we set the local coordinate system such that the z axis 
coincides with the equilibrium magnetization direction, 
M°' x = 0,M°' V = 0,M°' Z = 1, these correlations disap- 
pear and we have the same thermal fluctuations in x and 
y directions but no fluctuations in z direction, 
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and 



= 0. 



(29) 



Thus, the torque fluctuations produce effectively correla- 
tions and different values of thermal fluctuations in all 
other systems of coordinates different from the global 
one, where one of the axes is parallel to the equilibrium 
magnetization direction and where the equation of mo- 
tion for this component disappears. 

It is customary to introduce thermal fluctuations in 
the field components (see Q and originally W. F. Brown 
H) instead of the torque fluctuations as derived above. 
This has its origin in the representation of the LLG equa- 
tion in a spherical system of coordinates in form of Eq. ^. 
However, in both of these papers above only the axially 
symmetric case without interactions was considered. The 
big difference between these two approaches is the mul- 
tiplicative character of the field noise versus the additive 
noise of the torque. This turns out to be important for 
larger magnetization deviations. But first we will show 
that in the global coordinate system both approaches, 
torque and field, give the same result, as long as the 
magnetization deviations from the equilibrium are small. 

Let us use a decomposition of the field components 
according to H 1 — > H 1 + £ 1 , where are the components 
of the fluctuation part of the field. When this is done 
we obtain the following expansion of the equations of 
motion, 



~dV 



-e i3k M j H k - aH m [M m M l - 5 mi ] 

-e ijk M j £ k - of [M m M l - 8 mi ] 
A\M n ,H l ) + B ij (M n )^. (30) 



Furthermore, in the global system of coordinates we lin- 
earize the magnetization by the decomposition M % — > 



Ml, 



where 



are small fluctuations around the 



equilibrium values Mg, and apply the constraint condi- 
tion, |]Vi| = 1. For simplicity below we drop in the formu- 
las the particle index i. The components in the specified 
coordinate system are then 



dm x 
dm y 

~~dT 

dm z 



A x {rti)-{m y +am x )C + f x , (31) 
A y {rft) + (m x ~am y )C +f y , (32) 
K+m^f-O, (33) 
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where A (rft) stands for the linearized deterministic part 
of the LLG equation and 



(34) 



The constraint condition implies that in a first order ap- 
proximation it is m z (r) = 0, V r. This is compatible with 
Eq. |33| only if the field fluctuations can be considered 
to be small quantities, in which case products of the £ z 
with the m l can be ignored in Eq. (pT|) . These equations 
suggest that the field fluctuations contribute additively. 
From Eqs. 34 one can also obtain Brown's formulas for 
the field fluctuations (Eq. ||). 

It is important to note that the last equation is satisfied 
for any fluctuation field value due to the character of the 
LLG equation. Thus the f z value (or £ z ) is in this case 
undefined. In any case the component £ z is not efficient 
since it acts parallel to the magnetization direction. The 
assumption made in the paper of A. Lyberatos and R. 
Chantrell Q is that the field components are isotropic 
and that 



(35) 



This assumption in the global system of coordinates 
(where the fluctuation-dissipation theorem is applied) 
leads to the remarkable symmetry (|3^) of the field com- 
ponents in all the systems of coordinates and to the ab- 
sence of correlations. Furthermore, it is assumed that 
this property is valid through the magnetization rever- 
sal. 

For the torque fluctuations the reasonable hypothesis 
to mimic the field ones would be the assumption that 
there are never torque (force) fluctuations along the mag- 
netization direction. In this case the correlations between 
different noise components would appear in all other sys- 
tems of coordinates different from the global one. While 
equivalent near the equilibrium, these two approaches 
will be different far from it. At this point, we would like 
to restate that the whole theory is valid for small fluc- 
tuations around the equilibrium where both approaches 
coincide. 

In conclusion, the application of the Brownian dynam- 
ics approach to the motion of a magnetic system shows 
that interactions do not introduce correlations into ther- 
mal fluctuations introduced as both, either a fluctuating 
torque or a fluctuating field. Correlations may appear 
between different magnetization components as a result 
of the conservation of the value of the magnetic moment. 
The reasonable hypothesis that all the fluctuating field 



components are equivalent leads to Brown's well-known 
formulas for the fluctuating fields values without correla- 
tions. This validates all previously done micromagnetic 
calculations where this kind of assumption was made. 
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